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ABSTRACT 

Using optical-optical and optical-NIR colors, we analyze the radial dependence of age and metal- 
licity inside massive (M* > 10 10 - 5 M Q ), low-redshift (z < 0.1), early-type galaxies (ETGs), residing 
in both high-density group regions and the field. On average, internal color gradients of ETGs are 
mainly driven by metallicity, consistent with previous studies. However, we find that group galaxies 
feature positive age gradients, V*, i.e. a younger stellar population in the galaxy center, and steeper 
metallicity gradients, compared to the field sample, whose Vt ranges from negative in lower mass 
galaxies, to positive gradients at higher mass. These dependencies yield new constraints to models 
of galaxy formation and evolution. We speculate that age and metallicity gradients of group ETGs 
result from (either gas- rich or minor-dry) mergers and/or cold-gas accretion, while field ETGs exhibit 
the characteristic flatter gradients expected from younger, more metal-rich, stars formed inside-out 
by later gas-cooling. 

Subject headings: galaxies: clusters: general — galaxies: evolution — galaxies: fundamental parameters 
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1. INTRODUCTION 

Being hot dynamical systems, ETGs are often assumed 
to be formed from major mergers, where the progenitors 
can be either gas-rich disks or ETGs. The predominantly 
old stellar populations found in these galaxies indicate a 
rapid and intense star formation history, although the 
detection of comp act, massive ETGs at high redshift 
(|Daddi et al.l 120051 ) implies a chan ge in size of a facto r 
3-5 between z~2 and z=0 (see e.g. iTruiillo et alj f2007). 

Recent studies of the size evolution of massive ETGs 
over the past 8 Gyr propose a nu mber of alterna- 
tive scenarios, such as minor mergers (|Naab et al.l[20ol 
iTrujillo et alJl201l"i ) or a pufnng-u p of the central regions 
caused by the b aryonic mass loss (|Damianov et al.ll2009t 
iFan et aLll2010D . Alternatively, cold accretion - perhaps 
remin iscent of the monolithic collapse scenario (jLarsonl 
119751 ) - may provide an independent channel for the for- 
mation of massive galaxies, accounting for the highly ef- 
ficient st ar formation mea sured in these systems at high 
redshift (Dckcl et al. 2009). In addition, the environment 
can contribute as well to the star formation history. As 
a galaxy enters a group, its hot gas reservoir can be re- 
moved, shutting off star formation ("strangulation"). All 
these processes may be discerned if spatial information 
within galaxies is taken into account. The majority of 
studies on the star formation history of ETGs use obser- 
vations integrated within an aperture, losing the discrim- 
inating power. Radial gradients of photo-spectroscopic 
data represent a fundamental observable to disentangle 
the different formation channels described above. Studies 



of local samples of ETGs show that the majority feature 
red cores (see e.g. iPeletier et al.lll990D . Blue cores are 
less fre quent, being mo stly low-mass (M+<10. 105 Mq) 
ETGs ()Suh et all 120101 ) . Radial gradients of spectral 
line strengt hs reveal a significant trend in metallic- 
ity (see e.g. iForbeSj Sanchez-Blazquez. fc Proctorl 120051: 
Qgando et all 120051: iSpolaor et al.l 120091: iKoleva et all 



20111 ) and an intri guing dependence of the age gradients 



on e nvironment (Sanchcz-Blaz auez. Gorgas. fc Cardiell 
120061 ) . Extending the analysis to moderate redshift takes 
advantage of the lookback time, as younger populations 
are less affected by the age-metallicity degeneracy. The 
analysis of the color gradients of the GOODS sample of 
ETGs over a wide redshift range (0.4<z<l) confirmed a 
clear metallicity radial gradient for the red cored galax- 
ies, and a significantly you nger population at th e centers 
of the blue-cored galaxies ([Ferreras et al.|[2009h . with an 
increase in the fraction of blue- vs. red-cored profiles 
towards lower stellar masses. 

Galaxy mass and environment are the two main 
drivers of the format ion history of galaxies (see e.g. 
IWeinmann et all I2009D and recent work on ETGs has 
revealed a significant difference between the two, with 
environment onl y playing a second ary role on the stel- 
lar populations ()Rogers et alll2010l) although with more 
significant effects on dynamical relationships s uch as the 
Fundamental Plane (|La Barbera et all I2010al hereafter 
LLD10) . In this letter, we extend the analysis by infer- 
ring age and metallicity radial gradients of low-redshift 
ETGs from a multi- wavelength dataset in the optical and 
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NIR. Our findings reveal a significant contribution of the 
environment to the formation and assembly mechanisms 
of ETGs. 

2. SAMPLE 

The SPIDER survey consists of a volume-limited sam- 
ple of 39,993 bright ETGs (M r < -20) in the redshift 
range 0.05 < z < 0.095, with griz photometry from 
SDSS-DR6. 5, 080 of these galaxies also have YJHK 
imaging from UKIDSS-Large Area Survey ( see LLDI0, 
and references therein, for details on data analysis and 
sample selection). ETGs are defined as galaxies with a 
prominent bulge], and a passive spectrum in their cen- 
tral regions (within the SDSS fiber). Galaxies have struc- 
tural parameters, i.e. effective radius, i? e , mean surface 
brightness within that radius, (/i) e , and Sersic index, n, 
homogeneously measured from g through K with the 
software 2DPHOT, by fitting galaxy images with seeing- 
convolved 2D Sersic models. The environment of ETGs 
is characterized by an FoF catalogue of 8, 083 groups. 
A s hifting gapper techn ique is applied to this catalogue 
(see lLopes et al.l l2009al ) allowing galaxies to be classi- 
fied as either group members (~ 46%), non-group mem- 
bers (hereafter "field" galaxies; ~ 33%), or unclassified 
(~ 21%; see LLD10). Each group galaxy has 2D local 
density, Sjv, estimated by its distance to the n-th near- 
est group member, where n scales as the square root of 
group richness. 

We select only galaxies with better quality structural 
parameter s from q through K (~ 90% of the sam- 
ple; see lLa Barbera et al.l l2010bl hereafter LDDI0), 
and stellar mass M* > 3 x 10 10 Af©, which gives an 
M*-complete sample of ETGfl These selections result 
into a sample of 9,285 (936) field galaxies with opti- 
cal (optical+NIR) data available. For group galaxies, 
we select only ETGs residing in the highest density re- 
gions, Ejv > 16.6gals/Mpc 2 (E N > 6.4gals/Mpc 2 ) for 
the optical (optical+NIR) sample. These thresh- 
olds are chosen to emphasize the environmental depen- 
dence of galaxy color gradients, without overly reduc- 
ing the number of group galaxiet0, which consist of 
1,792 (optical) and 487 (optical+NIR) ETGs. Most 
of the parent systems where these ETGs reside are rel- 
atively poor groups, with an average velocity dispersion 
of ~ 370 ferns -1 . Both field and group ETGs are split 
among three M*-bins: "low" (3<M*<5xl0 10 M Q ), "inter- 
mediate" (5<A/*<8xl0 10 M Q ), and "high" (8<M*<45x 
10 10 M Q ), being the number of group ETGs the same in 
each bin. 



3. INTERNAL COLORS VS. GALAXY MASS AND 
ENVIRONMENT 

The internal color gradient of an ETG, defined as the 
slope of its radial color profile, written as X7 g -x, where 
X = rizYJHK, is estimated as described in LDD10, 

1 accounting for > 80% of total light, 

2 The Mi, limit is established with the same approach as in fig. 3 
of LDD10. 

3 Since optical-NIR color gradients do not depend significantly 
on Ejv (Sec. |3j, using the same Sjy threshold for both optical 
and optical+NIR samples would not change our results. However, 
adopting a lower Ejv threshold for the optical+NIR sample allows 
us to effectively reduce error bars. 
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Fig. 1. — Color gradients of ETGs as a function of environ- 
ment and stellar mass. Group and field ETGs are plotted with 
empty and filled symbols, respectively. Low-, intermediate-, and 
high-mass galaxies are plotted with squares, triangles, and circles, 
respectively. Upper panel: median V 9 _ r versus median stellar 
mass. Error bars are 1 <r confidence intervals on median values. 
The inset shows how central values of g — r, (g — r)o.i, depends 
on stellar mass and environment. Lower panel: the same as upper 
panel but for the optical-NIR color g — K. 



using deconvolved Sersic models of each galaxy to mea- 
sure g — X colors on concentric ellipses with axis ratio 
and position angle fixed to the r-band Sersic fit. A line is 
fitted to the resulting g — X vs. lo g R profile, in the ra - 
dial range of 0.1 to lR e (as in, e.g. iPeletier et alJ l990). 
where R is the equivalent galactocentric distance and R e 
is the r-band effective radius. The line slope gives V g _x- 
Fig. [T] shows the dependence of median optical-optical 
and optical-NIR color gradients, V g _ r and V ff _K0, as 
well as central colors, (g — r) .i and (g — K)o.i (com- 
puted at R = 0.1i? e ), on galaxy mass and environment 
(i.e. field vs. group). To a first approximation, optical- 
optical colors are more sensitive to the age of a stel- 
lar population than its metallicity, while optical-NIR 
colo rs are more s ensitive to metallicity than age (see 
e.g. Ide Jond ["19961. Hence, the upper and lower pan- 
els of Fig. [1] can be roughly seen as reflecting the be- 

4 The other available colors (g — i, g — z, g — Y , g — J , and 
g — H) are not shown for brevity, but are included in the stellar 
population analysis of Sec. [4] 
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havior of age and metallicity, respectively. At a given 
stellar mass, group galaxies have significantly shallower 
optical-optical color gra dients than their fi eld co unter- 
part^l, consistent with lLa Barbera et all j|2(ll)") i. who 
found cluster galaxies to have systematically flatter V s _ r 
(by ~ 0.03mag/dex at z < O.f) with respect to galaxies 
in less dense environments. However, optical-NIR color 
gradients do not exhibit any significant environmental 
dependence, implying that age (rather than metallicity) 
gradients are changing with the environment. The small 
panels in Fig.[T]also show that, at all masses, group galax- 
ies have redder central colors than those in the field. On 
average, the color difference amounts to ~ 0.02 mag in 
g — r, and ~ 0.05 mag in g — K. This may result from 
both age and metallicity increasing with mass, and group 
(relative to field) ETGs having ol d er ste ll ar populations, 
as found by, e.g.. iThomas et al.l ([2005D ; IGallazzi et al.l 
(poMk iPasauali et al.l (|2010D . 

4. AGE AND METALLICITY GRADIENTS VS. GALAXY 
MASS AND ENVIRONMENT 

We assume that the dependence of internal col- 
ors of ETGs on mass and environment is mainly 
driven by the age and metallicity of the underly- 
ing stellar populations, which is motivated by the 
existence of well-establi shed absorption-line gr a dients 
in ETGs (see, e.g., iSanchez-Blazquez et al.l 120071: 
iRawle. Smith, fc Lucevl I20f0f ). Although a dust com- 
ponent may also generate color gradients ()Silva fc Wise! 
f996), no observational eviden ce has been found so far 
that dust can play a major role (jSavov. Welch. Fich]|2009t 
LDDfO). 

We fit simultaneously all available central colors and 
color gradients (as shown in Fig. [IJ with stellar popula- 
tion (SP) models, inferring galaxy central age and metal- 
licit j0, as well as age and metallicity g radients Q, Vt and 
V z (as in, e.g., iTortora et all 12010). Each SP model 
cons ists of a p air of Simple Stellar Populations (SSPs) 
with lChabrier) (|2003| ) IMF, from the Chariot & Bruzual 
(201 1, in preparation; CB11) synthesis code. The "in- 
ner" SSP is used to model central galaxy colors, while 
the difference of g — X between the two SSPs is fitted to 
V 9 _X) providing Vt and V z- The fitting is done in a \ 2 
sense by varying age and metallicity of both SSPs in the 
range 1 < t < 14 Gyr and 0.2 < Z/Z Q < 2, respectively. 
We found that simila r resu lts are obtained when us- 
ing |Bruzuar^£harlo3 (|2003f ) population synthesis mod- 
els, or using exponentially declining star-formation mod- 
els. As shown in Fig. ® the age and metallicity content of 
ETGs depends significantly on both galaxy mass and en- 
vironment, reflecting the behavior of optical-optical and 
optical-NIR colors of Fig. [T] For all environments and 
mass bins, the main driver of color gradients in ETGs 
is a significantly neg ative (< —0-3) metallicity gradient 
(consistent with e.g. iPeletier et al.| [T990; Fcr reras et al.l 

5 We note that this trend cannot be caused by an increase of 
the E vs. SO fraction in the gr oup environment, as l enticulars have 
shallower color gradients (e.g. [Rocdigcr ct al. 2 0111) . 

6 Central age and metallicity are not discussed here, as their 
estimates are less robust than that of SP gradients, involving an 
absolute, rather than a relative, matching of model to observed 
colors. 

7 i.e. The variation of logarithmic age and metallicity per decade 
in galactoccntric distance. 
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Fig. 2. — Age (upper-panel) and metallicity (lower-panel) gradi- 
ents of ETGs as a function of stellar mass and environment. Sym- 
bols are the same as in Fig. [l] Error bars, defining 1 a confidence 
intervals, are obtained by shifting the (g — A")o.i's and V s _x's 
according to their errors, and repeating the fit. 



2009) , with field galaxies having somewhat shallower V z 
(by ~ 0.06) than their group counterparts. However, age 
is found to play a significant role, and, perhaps more im- 
portantly, this depends on the environment where galax- 
ies reside. Group galaxies have positive, albeit small , 
age gradients (~ 0.1, i.e. an age variation of ~ 23% 
per radial decade), with the galaxy stellar populations 
being younger in the center than in the outskirts, while 
field ETGs exhibit systematically smaller, mostly nega - 
tive, age gradients, consistent with iFerreras et all (2009) 
based on the analysis of the evolution of the color gradi- 
ents of field ETGs from GOODS data. Notice, however, 
that at high-mass field galaxies also tend to have slightly 
positive V t (~ 0.04 ±0.01). 

5. DISCUSSION 

The analysis presented in this Letter suggests that 
there is a significant trend of the internal distribution of 
age and metallicity in massive ETGs with galaxy mass 
and environment. The results of Fig. [2] can be qualita- 
tively discussed in the light of current formation scenarios 
of ETGs, considering that most (~ 90%) of our sample 
of group galaxies consists of satellite, r ather than centra l, 
galaxies (according to the definition of lYang et al^ OO?): 
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- Gas-rich mergers. The presence of a younger stellar 
population (i.e. a positive age gradient) in the center 
of group ETGs supports a dissipative formation picture, 
whereby gas-rich mergers fuel the central region with 
cold gas, allowing for the formation of younger, metal- 
enriched, stars. In a hierarchical picture of galaxy forma- 
tion, the more massive ETGs start formi ng stars at ear- 
lier e pochs (albeit assembling later, see iDe Lucia et al.l 
2006), when disk-like pro genitors likely had a mo re tur- 
bulent and clumpy ISM (|Forster Schreiber et aLll2009f ). 
possibly implying a larger amount of dissipation. This 
might explain why age (metallicity) gradients tend to 
marginally increase (decrease) with stellar mass. How- 
ever, this is not supported by SPH simulations of galaxy 
formation, with merging prod ucing no correl ation of 
metallicity gradients and mass (|Kobavashill2004| ) . More- 
over, both the light profile shape as well as the size- 
mass relation of ETGs at z ~ provide evidence 
against an increasing importance of gas dissipatio n at 
high mass (jHopkins et al.ll2009HShankar et al.H201lD . 

- Cold accretion. An alternative scenario to explain 
positive age gradients involves the accretion of cold 
(clumpy) gas from the cosmological surroundings, lead- 
ing to subsequent star formation in the central region. 
Cold accretion is expected to be more important at high 
rcdshift (z > 2), being a relevant mode for the forma- 
tion of galaxies as massive as ~ 1O 12 M (IDekel et al.l 
2009, but see also Ker es et ail I2005D . In a cold ac- 
cretion scenario, where the formation process resem- 
bles monolithic collapse, the Vt and V ' z would natu- 
rally correlate with galaxy mass as star formation lasts 
longer in the center of more massiv e systems having 
a deeper central p otential well (e.g. IKobavashl I2004t 
iPipino et al.l l2010f ). This interpretation is also sup- 
ported by the fact that: (1) the ages and metallicities 
of the stellar populations within and amon g ETGs are 
tight ly correlated to the local escape velocity (jScott et al.l 
2009); and (2) the resulting gradients depend on mass 
(e.g. iForbes. Sanchez-Blazquez. k, Proctorll2005D and, at 
fixed mass, on the duration of star- formation, a property 
that can be parametrized by the a to iron abundance ra- 
tio ( LDD10). One should notice, however, that the ex- 
istence of a correlation between Vz and mass has been 
longly debated in the lit erature (see, e.g., ISpolaor et al.l 
2009; iKoleva et aT]|201lL and references therein). 

- Later gas accretion. While group galaxies have 
their star formation quenched when entering bigger halos 
(through "strangulation" ) , field galaxies can accrete gas 
longer, by radiative cooling of their dark-matter hot-gas 
reservoir, and form a younger, more metal-rich, stellar 
component outwards. This would lower, and eventually 



invert, the age gradient, making also metallicity gradi- 
ents shallower, consistent with what we see for field ETGs 
in Fig. H 

- Stripping. Environmental effects, such as ram 
pressure and/or tidal stripping, may also trigger (cen- 
tral) star formation in rec ently accreted group satellites 
(e.g. iBekki fc Couchl 120031 ). In general, stripping should 
be less important for high-mass galaxies, inconsistent 
with the trends seen in Fig. [2] 

- Dry merging. ETGs can also form by gas-poor in- 
teractions. Minor dry mergers and stellar accretion can 
drive the formation of the outer envelopes o f ETGs, ex- 
plain ing their size evolution with redshift (|Naab et al.l 
2009). Gas -poor mergers should mix the stellar popu- 
lations within galaxies , flattening p re-existing metallic- 
ity and age gradients (White 1980). Since V t and Vz 
steepen with mass (see also LDD10), our data seem to 
reject major dry-mergers as the main channel for the 
format ion of ETGs, unless the amount of flattening is 
small (jHopkins et al.ll2009D . Minor dry -mergers would 
lead to a deposit of metal-poor, old, stellar material in 
the outer regions of a galaxy, increasing (decreasing) 
the age (metallicity) gradients. Therefore, the trends 
of Vt and Vz in Fig. [2] might also result from minor 
dry- mergers to be more important at high-mass (con- 
sistent with model predictions; see iHopkins et al|[2010l; 
IDe Lucia et al.ll20Tlj ) and preferentially in group, relative 
to field, environments. However, such a trend would be 
at odds with the finding of a lack of environmental depen- 
dence of the redshift evolu tion of ETGs on the mass-size 
plane (|Rettura et al.ll2010[) . 

In summary, although we find that several mechanisms 
can contribute to the observed age and metallicity gradi- 
ents in ETGs (among them, cold accretion and minor dry 
mergers are the most likely ones), we do indeed detect 
physical differences in the formation mechanisms with 
respect to environment, with field (relative to group) 
galaxies having the characteristic flatter gradients ex- 
pected from inside-out star formation induced by late 
gas-cooling. Studies of radial gradients outside of the ef- 
fective radius will shed more light into the mechanisms 
that drive the formation of early-type galaxies. 
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